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We present a formalism for the matter effects in the Earth on low energy neutrino fluxes which is 
both accurate and has all advantages of a full analytic treatment. The oscillation probabilities are 
calculated up to second order term in e{x) = 2V {x)E / /^m"^ where V{x) is the neutrino potential 
at position x. We show the absence of large undamped phases which makes the expansion in e well 
behaved. An improved expansion is presented in terms of the variation of V{x) around a suitable 
mean value which allows to treat energies up to those relevant for Supernova neutrinos. We discuss 
also the case of three-neutrino mixing. 
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I. INTRODUCTION 



The_propagation of low energy neutrinos in the Earth 
[iLp , Q is an important aspect of phy sics of solar Q - 
nil and Supernova (SN) neutrinos lUl - ll^. It will be 
useful in determining the oscillation parameters, and, in 
future, to search for effects of 1-3 mixin g |l4l and for a 
'tomography' of the Earth (see, e.g. |23,|2ll). It might 
even be possible to look for small structures of the density 
profile [23. 

In the existing calculations of Earth matter effects (see, 
e.g. P] - [231 ) the density profile is often approximated 
by one, two or several layers (mainly mantle and core) 
with constant densities or a direct numerical integration 
of the evolution equation is performed. However, the 
emergence of the large mixing MSW solution to the solar 
neutrino problem opens a more efficient approach to the 
oscillation effects in the Earth. Indeed, for the LMA pa- 
rameters, the oscillations of the solar and (lower energy) 
supernova neutrinos inside the Earth occur in a 'weak' 
regime, where the matter potential V is much smaller 
than the 'kinetic energy' of the neutrino system, i.e. 



V{x) < 



2E 



(1) 



In ref . |23l | , the e perturbation theory was formulated in 
the basis of neutrino mass states Vmass = {i^ii ^2)^ . The 
oscillation probabilities and the regeneration factor were 
calculated to first order in e. The expressions obtained 
are valid for arbitrary density profiles with sufficiently 
low density . They simplify the numerical calculations 
substantially and allow to understand in details all fea- 
tures of the oscillation effects. The method reproduced 
immediately the analytic result obtained in 24] for an ap- 
proximate but realistic density profile. Similar integral 
expression for the regeneration factor has been discussed 
in 25]. 

Since e(x) increases with energy, the lowest approxi- 
mation in e{x) may not be enough for larger energies. 
For instance, if ~ 50 MeV (possible for SN neutrinos), 
WG find e(x) ~ 0.6 at the center of the Earth. 

The purpose of this paper is to improve on this method 
and obtain accurate formulas which are valid for higher 
energies. In section 2 the oscillation probabilities are cal- 
culated in second order in e{x) and the convergence of 
the e expansion is commented on. In section 3 we sug- 
gest an improved perturbation theory which allows one 
to extend the expansion to higher energies. The general- 
ization to three neutrinos is given in section 4 and a brief 
conclusion in section 5. 



Here V{x) = \/2GFNe{x), Gp is the Fermi constant, 
Ne{x) is the number density of the electrons, Am^ = 
m\ — m\ is the mass squared difference, and E is the 
neutrino energy. 

In this case one can introduce a small parameter 



e(x) 



2EV{x) 
0.02- 



(2) 



E 




r^e(x)i 




10 MeV 









7.7- 10-5 eV^ 



II. SECOND ORDER CORRECTIONS TO THE 
OSCILLATION PROBABILITIES 

In this and the following section we consider the mixing 
of two (active) neutrinos Vf = U{9)i'mass, where Vf ^ 
{ve,Va)^ and Vmass = (1^1,1^2)'^ are the flavor and mass 
states, respectively and Va is a linear combination of i/^ 
and Vr- U{9) and 9 are the mixing matrix and mixing 
angle in vacuum. We define the matrix U{a) as 



where Na is the Avogadro number, and consider an ex- 
pansion of the oscillation probabilities in t{x). 



U{a) 



cos a sma 
- sin a cos a 



(3) 
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In 1231 the following expression for the S'-matrix in the 
mass eigenstates basis was derived p^ : 



S 







"V 


dx ^ 


J Xo 








— dx 


jdy 


Jxo 


J Xo 


+ ■•■ , 





1 

e^^"o^^ J + 



T(x) 



e 
1 



1 

e*"*"o- 



where 



dxA"(a;) 



(4) 



(5) 



is the adiabatic phase difference acquired by the neu- 
trino eigenstates in matter on their trajectory between 
two points xi and X2. A™ (a;) is defined as 



A'"(a;) 



2E 



Vl -2e(x)cos26l + e(x)2 ; (6) 



in vacuum we obviously have 



A" ^ A = 



Am^ 
'2E 



(7) 



The S'-matrix in Q is written as a perturbative expan- 
sion in T(a;) where 



T{x) 



sin 26* 



01\ 1 
loj + 2' 



V{x)[ — )+;lA™(x)sin2-'' 1 



(8) 



2 • V 1 / ■ 2~ V -f 

9' is the mixing angle of the mass eigenstates in matter, 
e sin 26* 



sin 26*' = 



(cos26'- e)2-t-sin2 29 



e sin 261™, (9) 



and 0™ = 6' -I- 0' is the corresponding mixing angle of the 
flavor states. 

The S'-matrix in eq. Q refers to a straight path through 
the earth from the entry point xq to an exit point Xf 
and the coordinate x is measured along the path. For 
notational convenience, we do not put labels xo,a;/, etc. 
on S. 

Using eq. ((HJ, we obtain the S matrix in terms of the 
potential V: 



S-- 



1 \ .sin 2^ rt^, J e"^"o- 



e 







sin^ 



'''' l)fZ.vi.f 



4A VO -e 

,2 



(10) 



sin2 20 r /e'^"-^^ 



The two last terms (proportional to e^) come from the 
first order in T (term proportional to sin^ 9' in Eq. 0) 
and the second order in T (see Eq. Q) correspondingly. 

Using the evolution matrix in the mass state basis I|1U|I . 
we can calculate the amplitudes and probabilities of var- 
ious transitions. The evolution matrix from the mass 
states to the flavor states relevant for the solar and SN 
neutrinos equals U S, where U is the vacuum mixing ma- 
trix (j3Jl. Consequently, the amplitude of the mass-to- 
flavor transition, is given by 



(11) 



The probability of the 1^2 ^ i^e transition, P^. 



I A. 



\U,,{9)S,2\ 



is then found to be 



sin^ 9 



sin^ 29 cos 29 



^sin^29 nxy(x)sinC-.., 



dx 



dyV{x)V{y) cos<_,(12) 



where the last term is the correction. The integra- 
tions over X and y can be disentangled. Indeed, writing 
(l^^-^x — 'f'y^-^z + "/"z^xj where z is an arbitrary point of 
the trajectory, we find 



dx f dyVix)V{y) coscf)'^ 

Xo Jxo 



m 



dx V{x) cos 



+ 



dx V{x) sin 0" 



(13) 
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This shows that the second order correction is positive 
for all V which do not vanish. 

Furthermore, for a symmetric density profile (with re- 
spect to the middle point of the trajectory) the second 
term in (|13() vanishes. This can be seen immediately by 
choosing z — x = (xf + xq)/2 in the center of the trajec- 
tory. So, finally we obtain for a symmetric profile 



= sm 



- sin^ 29 
2 



dxV{x) sin<?!)™_,^^ 



1 



■ sin^ 29 cos 29 



n 2 



dxV{x) cos 



(14) 



or (using again the symmetry of V) 



P. 



i sin2 26l sin (/.J 
1 9 

- sin^ 29 cos 29 



dxV{x) cos(/>" 



dxV"(a;)cos (/)™ 



(.15) 



The phase (t^^-^^f should be calculated according to 

The two last terms in (|15f) determine the regenera- 
tion parameter defined as freg = Pv2—>ve ~ sin^ 9 (see, 



3 



e.g., 23)- The probability of the Vi i/e oscillations 
can be obtained immediately from the unitarity condi- 
tion Pi/i^y^ = 1 - P^^^^^. 

According to (|15|l the effective expansion parameter of 
the series is 



dxV{x)coa 0™ 



so that 



p 

1/2- 



sin^fl + sin^ 20 



sm 



I 



cos 29r 



(16) 



(17) 



Notice that here the adiabatic phase should be calculated 
from the center of trajectory to a given point x, which 
corresponds to the explicit analytic expression obtained 
in Ref. (2^. According to (|17(l the first order correction 
is absent for trajectories with (t>™^xf = """^i = inte- 
ger) and the second order correction would be zero for 
maximal vacuum mixing. 

Taking A„i « A we obtain the useful bound 



2E 



dyVijj) cosy < 



2EVr,, 



Am^ 



(18) 



Vmax is the maximum value of the potential on the tra- 
jectory and y{x) — ^^x. 

In ea. ((Tn|l we note the presence of a possibly large 
phase (/)™_,^^ and an undamped integral in the term 
~ F(a:)^ (see 1-1 element of the matrix). It originates 
from term in T. (The undamped terms are absent 
in the linear term in e H^-) This could be a prob- 
lem, because the potential (squared) is integrated over a 
large distance without an oscillatory damping, and this 
might give rise to a large second order term in the ex- 
pansion. However by a simple partial integration of the 
last, ~ V{x)V{y), term in UlUI) one can see that the un- 
damped integral cancels. We have verified that this also 
happens in order for constant potentials. Therefore 
the e expansion appears to be well behaved (see also [26] ) . 



III. IMPROVED PERTURBATION THEORY 

As mentioned before, the accuracy of our expressions 
decreases for higher densities and energies. However, the 
expansion parameter can be reduced and therefore the 
expansion can be improved. This can be achieved by 
considering a perturbation around some average poten- 
tial Vb rather than around the vacuum value Vb = . 



In this case we expect the expansion parameter to be 



2El^V _ 2E(y - Vb) 
Am^ Am? 



(19) 



The corresponding results can be immediately obtained 
from the original perturbation theory. Indeed, the tran- 
sition to an average potential Vq is equivalent to consid- 
ering the problem in the basis v'^ = (j^i,J^2)i where v'j' 
are the eigenstates of the Hamiltonian in matter with 
a constant potential Vq. These states are analogous to 
mass eigenstates in the Vq — theory. Therefore the 5- 
matrix for (i/j^, 1^2) follows from the S matrix for mass 
eigenstates (|10|) by the substitution 



V ^ AV = V -Vq, 



(20) 



where 0™ is the flavor mixing angle in matter with the 
potential Vq: 



sin26'S 



sin 26* 



^l-2eo cos 20 + eg 



and 



£0 



2EVq 
Am^ 



(21) 



(22) 



The adiabatic phase differences generated for the eigen- 
states traveling in matter with true V are invariant un- 
der a shift of the average potential, so that the phases, 
4>T.^x- are unchanged. Therefore 



S° = S'(AF,6'^). 
We introduce the mixing matrix 
C^o = U{9'q) 



(23) 



(24) 



which relates the eigenstates of neutrinos in the potential 
Vo to the mass eigenstates in vacuum: Vmass — UqV^. 
The angle 9'q is given by 



sin 26*0 = eq sin 26*™ 



(25) 



and it is easy to check that 9 = 9"^ — 9q. 

Now the amplitude of the mass-to-flavor transition, 
i^i ^ i^a, equals 

A.,^.. =f/a,(C)(^").fct^L(^o)- (26) 
A straightforward calculation leads to the i>2 i^e 



oscillation probability P^, 
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P^,^^^ = sin^ 61 + eo sin^ 26'^" sin^ .1^1^ + gin^ 20J[' cos 29'^ / da; (x) sin (j)'^^^^ 

I sin' 29^ cos 200- rdxAV{x) sm^Z^, - '-^ sin^ 29^ Tdx /lyAy(x) Al^(y) cos(</>™_-</)-_,,^) 
1 r^j r^f 



- sin2 20j;^(cos26'J," + cos26'() - 2sin2 6'- eosin2 26'J,") dx dy AV (x) AV (y) cos c 

° J'XO Jxo 



>^L. • (27) 



r 



We note that there are two first order (in AV^) terms, 
one containing the other 4>™^^^ in contrast to the 

original theory which contains the phase (jj^-^x, only. 



For Vo = eq. H27|l coincides with the previous result 
For a symmetric density profile we obtain 



P^^^^^ = sin^ 61 + eo sin^ 20^ sin^ (^g^,,^^ + - sin^ 29"^ (cos 26*;, + eq cos 26*^) sin (^gl,^^ / dx AV{x) cos cfi^^^^ + 



+ - sin^ 20^ (cos 261™ + cos 29'^ - 2 sin^ 9 - 2eo sin^ 29'^') 



dxAV{x) cos0^!_^ 



(28) 



Thus the effective expansion parameter of the series in 
the improved perturbation theory is 

r'dxAV{x) cosC,. (29) 

J X 

The choice of Vq is arbitrary; the full expansion of the 
S matrix does not depend on it. It just should be chosen 
in a clever way. 

To illustrate the improvements, let us consider neutri- 
nos with energy 50 MeV [100 MeV]. For such neutrinos 
£=0.2 [0.4] in the upper mantle and e=0.6 [1.2] in the 
core. Thus, the average is eg — 0.4 [0.8]. Without im- 
provement, one expects the accuracy of the computation 
of ~ 0.2 [0(1)] in the core; with the improvement it 
is reduced to (e - eo)^ ~ 0.01 [0.06]. The optimal Vo 
can be chosen independently for each trajectory inside 
the Earth. For a mantle crossing trajectory, for instance, 
one would take the average value in the mantle. 

A 'good' value of eo may come from the observation 
that the second order term in H28|) is multiplied by the 
prefactor 



0.45 



0.4 



0.35 



0.25 



0.2 




0.2 0.3 0.4 0.5 0.6 



tan e 



FIG. 1: The dependence of eo = 2EVo/Am^ on tan^ S ob- 
tained by setting the prefactor in eq. (I3U1 equal to zero. 

IV. CORRECTIONS DUE TO 
THREE-NEUTRINO MIXING 



(cos 261™ cos 261^ - 2 sin^ 9 - 2eo sin^ 261™) . (30) 

Since eg is arbitrary, one may choose it such that this 
prefactor vanishes. In Fig.l we show the tan^ 61 depen- 
dence of eo for which the prefactor vanishes. 

In the limit V ^ Q the second, the third and the forth 
terms in (|28|l cancel each other ( up to ej^), and the prob- 
ability reduces to sin^ 9. 



In the standard parametrization the lepton mixing ma- 
trix is 

U = O23 diag{l, 1, e'^-)0i3 diag{l, 1, e~'^-)Oi2 = 

fCl3Cl2 C13S23 Si3e~**'^P 

-S12C23 - Ci2S23Sl3e^^'"' C12C23 - Si2S23Sl3e"'='' C13S23 
S12S23 - Ci2C23Sl3e"^<=f -C12S23 - Si2C23Sl3e'*'=f C13C23 , 
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(31) 



By a redefinition of the mixing matrix 
U ■diag{l,l,e'^'"') 
the Hamihonian becomes real, i.e. 

00 0\ /voo\ 
n ^ \ As +C/M000 C/ 
00 Aa/ \0G0/ 

Vc\^cl2 Vcl^Si2Ci2 VC12C13S13 

Vcl^Si2Ci2 As + VCI3SI2 ^Si2Ci3Si3 | ,(33) 
VC12C13S13 VS12C13S13 Aa + Vsjs 

where = Am%/2E and Aa = Aml^^/2E - As. 

Thus we see that both the CP phase Sep and 023 do 
not influence the propagation in matter (determined by 
the Hamiltonian). Also, since in (|31|l the first line does 
not contain Sep and 623 these parameters disappear in the 
oscillations from v^. to Ve, or from I'e to mass eigenstates 
and vice versa. They manifest themselves only when one 
considers the flavor states i^^ or i/^. 

These arguments are general and are valid in arbitrary 
matter density. 

We now write the Hamiltonian in the form 



where 



1-/ "T-ZO 



(31.), 



where 



(34) 



(35) 



and 

V+A +A —A"— A™ 
T(3,) = H-nl+diag{0, A„ A.) - I ' ° / 

sin 26*12/2 C12S13/C13 
= yc^3 I sin 2^12/2 S12S13/C13 I +0(F').(36) 




C12S13/C13 S12S13/C13 







A™ and A™ are the eigenvalues of the Hamiltonian in 
matter [2^. In eq. (|36|l we have subtracted a term pro- 
portional to the unit matrix in order to make it traceless 
and thus convenient for a power expansion. 

A straightforward calculation leads to the transition 
probabilities of the mass eigenstates to i^e- 



2 2 
'C13C12 



sin-" 261 



12 4 



Ci3 dxV sin (p 



-2 C22C]^3S]^3 / dxV ^IW 

'4^x — >-x f 1 



(37) 



2 2 



sin^ 26*12 



-13 



dx V sin ( 



r, 2 2 2 
^ ■'12C13S13 



dxV sin{'ipx->xr4'x->xf), (38) 



P.. 



SC-inC-ioSi 



'13 ^ ^ "-12"-13^13 / dxVsmil; 

X — >X f 

' xo 

+2si2C?3*i3 / dxYsmiipx^x^-cly^^j;^), (39) 

Jxo 



= AT{x) dx , V'a-.b = A^(x) dx . (40) 



a 



a 



(32) The function sin-tp^^^^ oscillates AJ^/A'^ 



Thus, 



Amltm/AniQ times faster than sine 
the corresponding integral is roughly Am^j^/Amg 
times smaller than the one which contains the phase 0; 
furthermore, it has a prefactor S13. Therefore we get to 
a good approximation 



P^,_,^ = 0^30^2 - —i7-^c^3 / dxVsm(t)^^^^,{'il) 



P 

-f^l/3^l/e— *13- 



2 2 '^i'^^ 26*12 4 



Ci3 / dxV sin(j)x^x^,{A2) 
(43) 



These results may be also obtained from eq. p5l | 
(see [13 for some earlier discussion). If Aa ^ As ^ V 
and S13 ^ 1, the third neutrino decouples and one arrives 
at the two neutrino propagation problem in matter with 
potential V — > ^Ci3 and mixing angle 6*12. Following the 
procedure of section II and using the full mixing matrix 
U diag{l, 1, e**'''' ) we easily recover eqs. (HTll - (|^ . 



V. CONCLUSION 

Motivated by the large mixing MSW solution to the 
solar neutrino we have developed a simple formulation of 
the earth matter effects on low energy neutrino beams. 
Following l23|, we derive an expansion for the neutrino 
transitions in terms of the parameter 



eix) 



2EV{x) 
Am? 



to second order. By choosing a convenient constant aver- 
age value for the neutrino potential as starting point, the 
precision can be substantially improved and it is possible 
to reach an accuracy of a few percents even for energies 
near 70 — 80 MeV. The effective expansion parameter is 
a simple integral in eq. (or eq. if^ ) together with 
eq. JSJ) which can be done numerically. The expansion 
allows for a convenient quantitative discussion of vari- 
ous physical effects such as the attenuation effect to the 
remote structures of the density profile or the effect of 
energy resolution of detectors. We also consider the case 
of three-neutrino mixing. 
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